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Abstract 

We investigate quantum tunneling in the theory of a complex scalar field with 
a global U(l) symmetry when the charge density of the initial configuration does 
not vanish. We discuss the possible final configurations and set up the Euclidean 
path integral formalism to find the bubble nucleation and to study the bubble 
evolution. For the stationary path, or the bounce solution, in the Euclidean time, 
the phase variable becomes pure imaginary so that the charge density remains 
real. We apply this formalism to examples when the initial charge density is small. 
While the phase transition considered here occurs in zero temperature, the bubble 
dynamics is richly complicated, involving conserved charge, the sound wave and 
the supersonic bubble wall. 
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1. Introduction 



Recently, there have been some interests in the first order phase transition in- 
volving nonzero global charge. The finite temperature effective potential and the 
phase structure in this type of the phase transition have been extensively stud- 
ied.' 11 However, how the phase transition proceeds has not been discussed in detail. 
Here we investigate the phase transition in a model which involves nonzero charge 
and a nontrivial bubble wall dynamics. This model is the theory of a complex 
scalar field with a global abelian symmetry. Even in the zero temperature phase 
transition, this model has a rather rich dynamics in the phase transition, depending 
on the initial configuration and the potential energy. Since there are only two field 
degrees of freedom in this model, this model can be rather easily approached ana- 
lytically and numerically. We hope that our toy model illuminates some aspects of 
the phase transition involving global charge and that some of insights gained would 
be applicable to the QCD phase transition and the electroweak phase transition. 

The general formalism of the Euclidean phase intergal involving nonzero charge 

[2] 

in our model has been developed sometime ago. While this formalism has been 
applied to wormhole physics, there has been no direct attempt to apply for the first 
order phase transition. We extend this formalism to the first order phase phase 
transition at zero temperature, following the standard formalism. One interesting 
aspect of our formalism is that the stationary path of the angle variable of the 
complex scalar field becomes pure imaginary. From this Euclidean path integral, 
one can find the bounce solution and calculate the bubble nucleation rate. In 
addition, one can gain some insights on the bubble evolution. 

The initial configuration we are interested in here is a homogeneous configura- 
tion which is classically stable but not quantum mechanically. The charge density 
of the initial configuration is nonzero and uniform. The final configuration after 
the phase transition, it turns out, could be more complicated than the configu- 
ration of the lowest potential energy. In our model, there could be an attractive 
force between charges and charges clump together forming Q-balls. [41 Thus, the fi- 
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nal configuration could be inhomogeneous with Q-balls floating in the symmetric 
phase. 

The Euclidean path integral allows us to calculate the imaginary part of the 
energy for the metastable initial configuration by the semiclassical method. The 
contribution to the path integral is dominated by bounce solutions. Contrasted 
to the usual case, the initial charge density now breaks the 0(4) symmetry to the 
0(3) symmetry. In this paper we focus the case of small initial charge density, 
where the bounce solution is close the 0(4) bounce solution of zero charge. When 
the charge density is small, we can look at the perturbative correction to the 0(4) 
symmetric solution. The current conservation equation in this background turns 
out to be a boundary value problem in the classical electrodynamics and can be 
solved in the thin wall limit. This leads to some insights on the current flow in 
the bounce solution and the deformation of the 0(4) symmetric bounce solution. 
This in turn leads to an understanding of the bubble evolution via the analytic 
continuation. 

When there is nonzero charge density and the initial configuration is metastable, 
there is always the sound wave of the speed less than the speed of light. When 
the bubble of a "true vacuum" is nucleated, it will expand. The bubble wall speed 
could reach the sound speed of the initial configuration in finite time, becoming 
supersonic. In addition, some bubbles speed up to the speed of light in finite time, 
which implies some sort of a new instability. 

The plan of this paper is as follows. In Sec. 2, we introduce the theory of a 
complex scalar field. We study the stability condition of the possible initial config- 
urations and discuss the final configurations we expect after the phase transition. 
In Sec. 3. we study the Euclidean path integral. We offer a formal discussion of the 
bubble nucleation and evolution. In Sec. 4, we examine in detail various examples 
when the initial charge density is small, and present the qualitative pictures of the 
bubble nucleation and evolution. In Sec. 5, we conclude with some remarks and 
questions. 
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2. Model 



We study the theory of a complex scalar field = fe/y/2 with a global U(l) 
symmetry. The lagrangian is given by 

C = \d^\ 2 - U{y/2<j>) 

= \w) 2 + \f 2 W?-u(f) (2,1) 

The global symmetry arising from a constant shift of 9 leads to the conserved 
current 

J M = -i(0%0 - <9 M 0*0) = f%9 (2.2) 

The total charge is Q — J d 3 xf 2 9. While there could be global strings in this 
theory, they seem not to play any essential role in our discussion and will be 
neglected here. 

In this paper we are interested in quantum evolution of a metastable initial 
configuration with nonzero charge density. To start, we should have a proper 
description of possible initial configurations. While some of them could be in- 
homogeneous and localized, here we will focus on the static homogeneous initial 
configurations (/, 9 = wt) with a given charge density p = f 2 9 = f 2 w. The energy 
density for these configurations is given by 

e(pJ) = ^ + U(f) (2.3) 

The first term in the right hand side of Eq.(2.3) is the centrifugal term due to 
the conserved charge. An initial configuration would lie at the local minimum of 
e(p, f) to be stable under homogeneous fluctuations, and satisfies 



| = -js + U'(f) = (2.4) 
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and Jjp > or 

o 2 

3^3 + U" > (2.5) 

We denote a local minimum by /o an d wo = p/fo- Fig.l shows an example of 
e(p, /) and C/ (/) for nonzero p. 

A metastable initial configuration should be stable, even under small inhomo- 
geneous fluctuations, fo + 5f and 9 = wot + 59. From the classical field equation, 
we get the linear equation for the fluctuations, 



- 5f + df5f + w 2 5f + 2w foS9 - U"(f )6f = 
-2w o f 6f-fi59 + f 2 d 2 59 = 



(2.6) 



With 5f,59 ~ e l ( at fe '^, Eq.(2.6) yields a dispersion relation 

[a 2 - k 2 } [a 2 -k 2 + wl- U"(f )} - 4a 2 wl = (2.7) 



or 

2 r2 3wZ + U" , l\3^ + C/" 2 , 2T2l 1/2 

a 2 = A; 2 + ^ I _ ±|(— Q — 2 ) +4w 2 k 2 j (2.8) 

One can see that there are one massive and one massless modes. The massless 
mode can be interpreted as the sound wave and satisfies the dispersion relation 
a 2 = v 2 k 2 for the long wave length, or for small k, where the sound speed is given 

vl = gm -f\ ( , 9) 



Thus, the configuration is stable under small fluctuations only if 

fSu"(f )-p 2 >0 (2.10) 
which is a stronger condition than Eq.(2.5). 
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When the charge density is high, the centrifugal term would be balanced by 
the highest order term in the potential, say U ~ ef n /n. From Eqs.(2.4) and (2.9) 
we see that the sound speed in this limit is given as 



Vs = J!LJ. (2.ii) 

V n+2 K ' 

It is interesting to note that the potentials of n = 4 in four dimensions, n = 6 
in three dimensions and n = oo in two dimensions are renormalizable interactions 
and that the corresponding sound speeds 1/a/3, 1 / y/2, 1 are those of hot relativistic 
gases in the corresponding dimensions. 

The stability condition (2.10) can be examined more concretely in the potential 

W) = ^/ 2 + f/ 4 + ^ 6 + ^ 8 (2.12) 

(Here we are not concerned about the renormalizability of the theory. We are 
interested in the general characteristics of the tunneling when global charge is 
involved.) Assume that the charge density is very small compare with other scale. 
When m 2 > 0, the condition (2.4) can be satisfied for /o near the symmetric phase, 

/0 = v / pM + O(p 2 ) (2.13) 
The stability condition (2.10) becomes 

foU'Vo) ~P 2 = 3s(-) 3 (l + 0(p)) > (2.14) 

TYli 

which is satisfied only if g > 0. The configuration is thus stable only if there is a 
short range repulsion due to the self interaction. The velocity of the sound in the 
case g > is given as 

^=B (i+o(p)) (2 - i5) 

which is much smaller than the unity when the charge density is small. 
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In the case where the potential (2.12) takes a local minima at / = v, U ~ 



which is automatically satisfied when the charge density is small. The sound ve- 
locity becomes 



This shows that the massless Goldstone boson becomes the sound wave as the 
nonzero charge density is introduced in the broken phase. 

We have examined the stability condition on the possible initial configurations. 
If an initial configuration is metastable, it will evolve quantum mechanically so that 
the potential energy is converted to the kinetic energy of bubbles and the radiation 
energy. The important question is how we know whether an initial configuration 
is metastable quantum mechanically? We would say a configuration is not stable 
quantum mechanically if we can find a field configuration of a lower energy and 
the same quantum number so that there is no superselection rule preventing the 
transition between these two configurations. 

In the case of the theory of a real scalar field, the answer comes immediately 
from the potential energy density. The potential energy density for the final con- 
figuration would be the lowest. A metastable initial configuration can decay via 
nucleations of bubbles whose interior is in the true vacuum. We know well how 
this phase transition proceeds. 

For the theory of a complex scalar field, the story is more complicated as 
there are two degrees of freedom which work together or against each other. If we 
regularize the system in a large finite box, we expect that the possible final config- 
uration has the minimum energy for a given total charge. The excessive energy of 
the initial configuration would be channeled into elementary excited modes in the 



^m ,2 (f — v) 2 with m' 2 = U"(v). The condition (2.4) implies p 2 ~ m /2 t> 3 (/o — v). 
The stability condition (2.10) becomes 



m'V + 3p 2 + vp 2 U"'(v)/m' 2 + C(p 3 ) > 



(2.16) 




(2.17) 
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final configuration through the radiations and the bubble collisions. In general, it 
is not easy to find such a final configuration. A possible final configuration could 
be inhomogeneous. In this paper, we concentrate on the cases where the potential 
U(f) has a local minimum at / = and the phase transition occurs between the 
/ Ri and / 7^ phases. Fig. 2 shows such potentials. In addition, we assume that 
the initial charge density is small compared with other scales in the problem. 

In these limits, we can analyze the problem in a somewhat satisfactory way. As 
far as the / field is concerned, this field wants to settle at the ground state of e(p, /) 
in Eq.(2.3). However, we will see that the charge will move so that the energy per 
unit charge e(p, f)/p takes the lowest value. Since we concern the phase transition 
between the symmetric and asymmetric phases, we choose U(f — 0) = 0. It turns 
out that there are three different cases of the phase transitions to consider: Case 
A) The ground state of U(f) is an asymmetric phase where U(f) < and the 
initial configuration is in the symmetric phase; Case B) The potential U (/) has a 
local minima at an asymmetric phase where U (/) > and the initial configuration 
is in the asymmetric phase; Case C) The potential U (/) is identical to Case B but 
the initial configuration is in the symmetric phase. In Fig. 2 U\(f) corresponds to 
Case A and U2(f) does to Cases B,C. 

Let us first study the first Cases B, C where the symmetric phase is the ground 
state of U(f). To find out the final configuration, let us recall that in these cases 
there could be Q-balls in the theory for an appropriate potential. Let us here 
recapitulate the Q-ball physics briefly.' 41 In the symmetric phase charged particles 
have mass m and so the ratio e/p would be m. For a given large charge Q, the 
lowest energy state however does not need to be made of a collection of these 
charged particles at rest, whose energy is m\Q\. Rather, that could be a Q-ball if 
there is an enough attraction between the charged particles. 

To find the condition on the potential for Q-balls to exist, let us examine 
a homogeneous configuration whose energy per unit charge is lowest. We first 
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minimize e/p with respect to p, 



dp 2/2 p- 
The charge density is then fixed as a function of /, 



d(e/P) - 1 -W> = 0. (2.18) 



p = 2/ 2 t/ (2.19) 



The energy per unit charge is given by 




(2.20) 



We have to minize e/p in Eq.(2.20) with respect to the / field. Fig. 3 shows \/2U / f 2 
for various potentials. When a/2£7 / f 2 takes the local minimum value w * at the 
nonzero f — f* field, such configuration is called Q-matters whose charge density 
is fixed to be p* = f%w*. A Q-ball is a sphere whose inside is made of Q-matters 
and outside is just the symmetric phase. When the total charge is large, the size 
of a Q-ball would be large and the surface energy will be negligible compare with 
the volume energy and the energy per charge would be very close to u>* Thus, if 
u>* is less than m, the Q-balls are stable against decaying into charged particles. 
Such Q-balls are possible with U\, Ui of Fig.3. 

For Case B, the initial configuration should be classically stable, which means 
g > as shown in Eq.(2.15). U2 in Fig.3 represents such a potential. 

When the stable Q-balls are possible, the final configuration for Cases B, C 
is inhomogeneous. There will be a region with Q-matters and the rest as the 
symmetric phase of zero charge density. The domain walls separating two regions 
would slowly evolve to reduce the surface energy. The ratio of the volumes between 
two regions would vary depending on the initial charge density. If the initial charge 
density is small, Q-balls will float in the symmetric vacua. If the initial charge 
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density is large, the balls of the symmetric vacua will float in the Q-matter. In 
Case B the phase transition would proceed with the nucleations of Q-balls and in 
Case C it would proceed with the nucleations of bubbles with less charge density. 

If Q-balls are impossible because w» > m as in Us of Fig. 3, then the final 
configuration would be in the symmetric phase with nonzero charge density. If 
the initial charge density is larger than p*, the final configuration could be more 
complicated would not be discussed here. 

To understand the condition (2.18) better, consider the pressure of a homoge- 
neous configuration given by 

P=jp-U{f) (2.21) 

We see Eq.(2.18) is identical to the zero pressure condition. When Q-balls can 
exist, we have argued that the final configuration could be inhomogeneous where 
Q-balls float in the symmetric vacuum. The zero pressure condition means that 
there is no pressure difference between Q-balls and the symmetric vacuum, leading 
to an equilibrium situation. 

Having analyzed the last two cases where the symmetric phase is the ground 
state of U (/), let us now consider Case A where the asymmetric phase is the ground 
state of U(f). For this case, both the / field and charge prefer the asymmetric 
phase because the the energy per charge in the asymmetric phase would be negative 
compare with that in the symmetric phase. The final configuration would be in 
the asymmetric phase with uniform charge density. An interesting observation on 
this case has been made when the initial state is not classically stable because 
g < in the potential (2.12). [51 Note that the initial configuration of zero charge 
density is classically stable. When we introduce a small charge density, charges get 
concentrated in several regions and these regions would evolve classically into the 
true vacuum. For the initial configuration at the symmetric phase to be stable, 
there should be a repulsive force between charges (g > 0) in long distance. 



10 



We note that the energy per unit charge and the pressure are sensitive to 
the shift of the potential energy. In the above arguments, U(0) = was used 
crucially because we are dealing with tunneling between the symmetric phase and 
the asymmetric phase. If we attempted to understand the tunneling between two 
asymmetric phases, we would have needed another device to figure out the final 
state. We will try to investigate this example elsewhere. 



3. Quantum Tunneling 

We have studied the general characteristics of the phase transition from a 
homogeneous initial configuration with nonzero charge density. We found out 
qualitatively what will be the final configuration after the quantum tunneling. 
We now want to approach this problem more analytically by using the Euclidean 
path integral formalism. When a nonzero global charge is involved, the standard 

[2] 

formalism should be extended to accommodate the nontrivial boundary term. Here 
we summerize and expand the known results. 

We start with the Euclidean generating functional, 

< F\e~ HT \I >= J [/d/d0]*£e _S|5 *j (3.1) 
where the Euclidean action is given by 

Se = J d 4 x |i(^/) 2 + \f\d,6f + [/(/)} (3.2) 
The initial and final states ^i t F describe the configurations of charge density, pi ; F, 



or 



*/,f(0, /) ~ exp ji J d 3 xp LF 9^ (3.3) 



In our case, the initial and final states are identical and describe the initial metastable 
configuration. 
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By calculating the imaginary correction to the energy of the metastable initial 
configuration, we can find out the tunneling rate or the bubble nucleation rate. 
When we sum over the multibounce contributions to the energy, 



e- £VT =< i\e- HT \i > 



K e- £oVT jl + VTKe- B + ^(VTKe- B f + .... j 
exp{-(£ - Ke~ B )VT) 



(3.4) 



where £q is the energy density of the initial configuration, Kq is the prefactor 
from the small fluctuations around the initial configuration, the B factor is the 
difference between the bounce action and the background action, SoVT, and KKq 
is the prefactor arising from the small fluctuations around the bounce solution. 
There is a negative mode around the bounce solution which implies the factor 
K pure imaginary. The bubble nucleation rate per unit volume is then given by 
2\K\e- B . 

The path integral will be dominated by the stationary configurations of the 
path integral. In our case the boundary condition on the / field of such stationary 
configuration is fixed to be the initial classical configuration. However, Eq.(3.3) 
implies that the boundary condition on the 6 field is free. Thus, the stationary 
configuration in the path integral (3.1) satisfies the Euler equation for the action 
Se + S, where £ is the boundary term 

£ = -i J d\ {p F 6(r, t f ) - Pl 9(f, 77)} (3.5) 

From the Euclidean equation from the combined action Se + £, one can see that 
the stationary path of the angle variable 6 should be pure imaginary, 

9 = — %t\ (3.6) 

The Euclidean field equations for the /, 77 field are then 

%f + f(d,v) 2 -U'if) = 

'V ./%'/) (l 
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The boundary condition on rj becomes 

f 2 ^(r,r I , F )=p I , F (r) (3.8) 

The / field should approach the time independent / given by the initial configu- 
ration. We assume that there is no vortex in the initial configuration. Then the 
classically stable configuration for given charge should satisfy 

%f+j5-U'(f)=0 (3.9) 

since d T f = and dii] at the boundary. In our case, the initial configuration is 
homogeneous in space and so Eq.(3.9) becomes identical to Eq.(2.4). 

The solution of Eq(3.7) is the so-called bounce for quantum tunneling. When 
the initial charge density vanishes, we know that the 0(4) symmetric solution of 
Eq.(3.7) has been shown to exist by the undershoot-overshoot method.' 31 In our case, 
the boundary condition reduces this 0(4) symmetry to the 0(3) symmetry because 
the charge density selects a preferred time direction. Thus, we are interested here 
the bounce which is 0(3) symmetric invariant under the spatial rotation. This 
makes Eq.(3.7) a partial differential equation, whose solution is much harder to 
find. We can either use some analytic tools or numerical analysis to find the 
bounce solution. In the next section, we use the perturbation method to get an 
approximate bounce solution when the initial charge density is very small. 

Once we find the 0(3) bounce solution we can calculate its action S F + <J 

from Eqs. (3.2) and (3.5). By using Eqs.(3.7) and (3.8), we can see the combined 
action becomes 

S E + <y = J d^x^d^) 2 + ^ft(d l r lh ) 2 + UUb)} (3.10) 
The bubble nucleation per unit volume is then given by Ke~ B where 

B = (Se + c) (bounce) — (S E + a) (background) (3. 11) 

While the action for the bounce Se + ct could be infinite, the difference B between 
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that of the bounce and that of the initial configuration should be finite when one 
expects a finite tunneling rate. While we will not attempt to calculate the effect 
of the fluctuations around the bounce solution, we note that the field fluctuations 
5f, 59 should be kept real in the functional integral. This is exactly what happens 
in a gaussian integral 



To find out the escape point or the bubble configuration at the nucleation 
moment, we use the time translation and reflection symmetries of the action under 
r — > — r and rj — > — rj of Se + v to choose the origin to be the center of the bounce 
so that 



As in our problem charge density remains real in Minkowski and Euclidean times, 
it is natural to identify the initial charge density of the bubble to be given by that 
of the bounce. The initial bubble configuration is then 



In usual quantum tunneling, momenta are imaginary and coordinates are real under 
the potential barrier, and we find the escape point in the coordinate space. In our 
case, the charge or momentum density and the angle variable have changed their 
role. Under the barrier the angle variable is imaginary and the charge density is 
real, and we find the escape point in the charge density. 

Once we know the initial bubble configuration (3.13), we can solve the field 
equation in the Minkowski time to find out how a given bubble evolve. The bounce 





(3.12) 



r/ b (r,T = 0) = 



f(r,t = 0) = f b (r,r = 0) 

e(f,t = o) = o 

%(r,t = 0) = 




(3.13) 
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solution (/ft, ?7ft) can also be analytically continued to a solution in Minkowski time, 



f(r,t)=f b (r,it) 

(3.14) 

9(r,t) = -ir) b (r,it) 

A further insight about bubble nucleation with global charge can be gained 

[2] 

by using the dual formulation. It is well known that Goldstone bosons can be 
described by an antisymmetric tensor field, B pv . In the Minkowski time, the dual 
lagrangian is given by 



Cdm = \{d,ff + ^jffjLp - U(f) (3.15) 

where H pvp = d p B vp + d v B pp + d p B^ v . The field strength of the antisymmetric 
tensor is related to the original current by 

f&B = l -e^H^ p (3.16) 

The uniform initial charge density becomes the condition of the uniform 'magnetic' 
field i?i23- In Euclidean time, there will be no boundary term arising from the wave 
function ty(f,Bij). The Euclidean lagrangian becomes 

Cde = \{d p f) 2 + ^plfyp + U(f) (3.17) 
The bounce equation becomes 

d lf + tL h U - u '(f) = 

D/ (3.18) 

We can relate the Euclidean fields either through the Euclidean time dual trans- 
formation or by comparing the current. The relation between the antisymmetric 
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tensor field and the angle variable in the Euclidean time is given by 

fd»r) = \e^°H vpa (3.19) 

Thus, the bounce solution in terms of the B^ v field would be real and there would 
be no contribution to the bounce action from the boundary. 

4. Examples 

We are now in position to examine in more detail three cases of the phase 
transition which we have discussed in Sec. 2. To be more specific we choose the 
potential to be given by Eq.(2.12). Case A has the initial configuration in the sym- 
metric phase. The potential energy has the absolute minimum at the asymmetric 
phase. For Cases B, C the potential energy U (/) has the absolute minimum at 
the symmetric phase and the local minimum at the asymmetric phase. The initial 
configuration of Case B is at the asymmetric phase and the initial configuration 
of Case C is at the symmetric phase. In Cases A, B the tunneling would proceed 
even when there is no initial charge density because the initial configurations are 
at the metastable points of U(f). Introducing a small amount of charge density 
would not change much of the original bubble nucleation. Thus we would expand 
perturbatively the bounce solution by the initial charge density and see how the 
zeroth order 0(4) symmetric solution deforms. These are the cases we will exam- 
ine closely in this section. For Case C, we do not have the zeroth order bounce 
solution because the quantum tunneling occurs solely due to the charge density. 
However, we can still get some insight for this case as we will see later. 

Let us consider first the cases where the initial configuration is unstable even 
without any charge density. The symmetric phase of the potential U\ and the 
asymmetric phase of the potential U2 in Fig. 2 are such initial configurations. When 
there is no initial charge density, the bounce solution can be obtained by the 0(4) 
symmetric ansatz, f(s = \fr 2 + r 2 )! 31 Let us assume that the thin wall approxi- 
mation works. We call that inside the wall, f = fi and outside the wall f — f e . 
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The wall radius a can be determined as follows. Suppose that the potential energy 
difference AU = U(f e ) — U(fi) > between f e and fi phase is small. Then the 
bubble radius will be large and we can approximate the bubble wall as a domain 
wall separating two phases. This wall satisfies the equation, d 2 j + U'(f) = 0, 
neglecting the potential energy difference. Define the tension of the wall to be the 
action density per unit three volume, T = J dx[(df / dx) 2 / 2 + U(f)}. The gain of 
the action due to this true vacuum bubble of the radius a is then 

7T 2 

S(a) = 2ir 2 Ta 3 - —AUa 4 (4.1) 

At a = 3T/AU, S(a) takes the maximum value 27tt 2 T 4 /2(A[/) 3 , which is the B 
factor in the bubble nucleation rate. 

We ask what happens to this 0(4) symmetric thin- wall bounce solution if we 
introduce a small initial charge density. From Eq.(3.7), we see that the equation of 
the angle variable is the first order in the charge density and the / field equation has 
a second order correction to the bounce equation of the zero charge density. Thus, 
we can solve Eq.(3.7) by a perturbative expansion around this 0(4) symmetric 
background. The phase variable will be first order in charge density and satisfies 
the current conservation 

di{P{s)d l r ] ) = (4.2) 

with the boundary condition f 2 r](f, r = ±oo) = po with the initial charge density 
Pq. The above equation can be interpreted as a boundary problem of a dielectric 
media in four space dimensions. The electric field is d^rj with the potential rj and 
the electric displacement is = f 2 d^r] with the dielectric constant f 2 . Eq.(4.2) 
implies the boundary condition at the wall that the normal component of and 
the tangential component of d^rj should be continuous. The boundary condition 
at infinity is that there is a constant external electric displacement field J T = pq. 
For a given 0(4) symmetric configuration described before, it is trivial to find the 
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potential rj e: i outside and inside the thin wall, 



Po ( 4/| 



Vi /e 2 V3/ e 2 + / 2 



(4-3) 



The charge density J T at the r = would be the charge profile of the bubble at the 
moment of nucleation as shown in Eq.(3.13). Let us now examine the implications 
of this solution (4.3) in various cases. The correction to the / field would be second 
order and will be considered in each case. 

Case A: from the symmetric phase to the asymmetric phase 

Let us first consider the case when the asymmetric phase is the ground state 
and the initial configuration is near the symmetric phase. Thus, f e ~ << fi ~ v . 
Since there is no initial charge density, the previous argument would imply that 
f e = 0. When we introduce uniform charge density in the initial configuration, 
the initial value of / would be given by Eq.(2.13) with 9 = mt, invalidating our 
assumption / ~ Oil). Here let us assume simply that / e is nonzero even when 
there is no charge, say due to a small bump in the potential at / = 0. This would 
not change the physics of tunneling under consideration much and allows us to use 
Eq.(4.3). 

The global current around the 0(4) symmetric bounce solution can be 
obtained from Eq.(4.3). Outside the thin wall (C > a), J^ — fld^e and inside the 
thin wall (£ < a), J M = f^d^rji. The energy per charge inside the wall is small and 
so the charge is attracted to the interior region, making the charge density inside 
the bubble be higher than that outside. The charge density profile of the bubble 
at the moment of nucleation would be given by J T at r = 0, 

o fl4 

"^exterior (1 ^4 )P0 ^ ^ 

"^interior 

when fi » f e . The charge density inside the bubble is four times larger than 
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the initial charge density. From Eq.(4.4), we can see J d 3 x(J° — po) = 0, implying 
that the charge inside the bubble came from the region near the bubble wall. Fig. 4 
shows this global current on the 0(4) symmetric solution background. 

Let us now consider the effect of the charge flow on the / field. The / field 
equation (3.7) can be expanded around the 0(4) symmetric solution as 

%Sf ~ U"(f)6f = - j§ (4.5) 

where J 2 can be obtained from Eq.(4.3) and / would be given by the 0(4) sym- 
metric solution. Rather than try to solve this partial differential equation, let us 
approach the problem more qualitatively. When f e << fa, one can show that at 
the north and south poles J 2 ~ 16pg both inside and outside and that at the equa- 
tor J 2 ~ outside and J 2 ~ f6pjj inside. The centrifugal term J 2 /2/ 2 would be 
then important at the poles but not at the equator. Directly from the / equation 
(3.7) and the previous argument about the thin- wall approximation, we see the 
centrifugal term reduces the tension on the domain wall by reducing the effective 
potential energy barrier. Since the tensions at the two poles are lower while the 
tensor at the equator remain constant, the curvature at poles would be larger than 
1/a and that at the equator will remain 1/a where a = 3T/AU is the radius of the 
0(4) symmetric shell. ( In addition AU in increased at the pole and remain un- 
changed at the equator, amplifying the curvature change between the poles and the 
equator.) Consequently, the 0(4) symmetric wall would be shrunk at the equator. 
We take a liberty to choose this wall configuration to be an ellipsoid 

r 2 r 2 

where b < a takes a complicated function of po- Fig. 4 shows the deformed bounce 
wall. 

Let us now think about the bubble nucleation and evolution. The bubble of 
true vacuum will nucleate with a radius b and a higher charge density 4po and then 
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expand. As argued in Eq.(3.14), we can analytically continue the bounce solution 
to find out how the bubble will evolve. From Eq.(4.6), we see that the bubble wall 
trajectory will be given by R(t) = y/|Vt 2 + ab. Since b < a, the terminal velocity 
vt = \Jbja is less than 1. 

How do we understand the finite terminal velocity? Energy conservation im- 
plies that the change of the wall energy comes from the potential energy difference, 

d (*jW) =d (** AUR3 \ (47) 



The tension of the bubble wall surface could depend on R. Integrating Eq.(4.7), 
we get 



R = dl-*£W (4.8) 
V AUR v ; 

When there is no charge density, T, AU are fixed and we see that the terminal 
velocity is the light speed. If the tension grows linearly with the radius for large 
R, T(R) ~ aR, the terminal velocity vt = \fbfa will be 



««=V 1 -|f <1 (49) 

The growth of tension, or the energy density of wall per unit be under- 

stood by considering the phase variable 9. While the charge density ffQ inside the 
bubble is larger than the charge density /| 9 outside, it is but not large enough to 
keep the phase variable space independent since fi » f e . The phase increases by 
mt outside the bubble and more slowly inside the bubble, leading to the increasing 
its space gradient at the bubble wall. This is what we think is the source of the 
increasing tension or energy density at the bubble wall. 

Let us now remind ourselves that the sound speed (2.15) in the symmetric 
phase is much smaller than the terminal speed (4.9) when the initial charge density 
is small. Thus, the bubble wall forms a sort of supersonic front in the symmetric 
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phase. However, our analysis is not accurate enough to compare the terminal speed 
and the sound speed (2.17) at the inside asymmetric phase. ( When the initial 
charge density is large enough, there is a possibility that the terminal speed is less 
than the sound speed in the symmetric phase. ) The thin wall approximation would 
fail eventually, because there is not enough charge lying outside the expanding 
bubble to keep the charge density inside the bubble to be four time larger than 
the initial value. The charge density profile around the expanding bubble should 
become more smoothly changing. 

Case B: from the asymmetric phase to the symmetric phase 

Let us now consider the case where the initial state is the broken phase which 
has the higher potential energy than that of the symmetric phase. When there is 
no initial charge density, there will be an 0(4) symmetric bounce solution, inside 
which the scalar field takes value near the symmetric phase, fi « f e . Again we 
ask what is the consequence of the small initial charge density. Eq.(4.3) implies 
how the current flows around this bounce solution. The charge density at the 
moment of bubble nucleation would be given by 

^exterior ~ U + g^JPO ^ 
interior ~ u 

when fi << f e . The charge is excluded from the symmetric vacuum region. Fig. 5 
shows the charge flow around the 0(4) symmetric bounce solution. 

The effect of the charge flow on the / field is given by the centrifugal term 
J^/2/ 2 as in Eq.(4.5). We can calculate J 2 for our bounce solution. Since f e >> fi, 
Eq.(4.3) implies that at north and south poles J 2 ~ both inside and outside 
and that at the equator J 2 ps inside and J 2 ~ outside. This raises the 
energy density of the false vacuum at the equator, lowering the barrier energy 
and increasing AU. This in turn lowers the tension of wall at the equator. The 
tension at the poles would remain unchanged and so the bounce solution would be 
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elongated at the equator. Now for the sake of the argument, we again approximate 
the wall as an ellipsoid, 

r 2 r 2 

^w- 1 < 4 - n > 

where b < a is a complicated function of p e . The deformed bounce solution is 
shown in Fig. 5. 



In the Minkowski time, the bubble wall trajectory is given by R(t) = ^a/b^/t 2 + b 2 , 
with the terminal velocity vt = \Ja/b > 1, which clearly violates the causality. 
Something should happens before the wall speed becomes the light speed. When 
there is a tachyonic mode, we say there is an unstable or exponentially growing 
mode. There are many possibilities. Since charge is pushed out from the bubble 
and is accumulated at the wall, the wall could stop expanding. Or the thin wall 
approximation could break down before the wall reachs the light speed. 

In Case B, the charge is pushed way from the bubbles and accumulated at the 
initial asymmetric phase. At the end of phase transition, we would be left with the 
islands of the original phase with high charge density, which are exactly Q-balls 
floating in the symmetric phase. 

Case C: forming Q-balls by quantum tunneling 

If the minimum of \f2U / f 2 is lower than the mass of the charged particles in 
the symmetric phase, charge likes to clump to Q-balls. Suppose the initial config- 
uration lies at the symmetric phase with very small charge density and is stable 
under local fluctuation. Since the minimum of the potential is chosen to be the 
symmetric phase, the initial configuration would be stable if there is no charge 
density. After small uniform charge density is introduced, the initial configuration 
however becomes unstable under the quantum mechanical tunneling transition to 
form Q-balls. Since we start from the minimum of the potential, we do not have 
the bounce solution at zero charge density. However, we can still gain some un- 
derstanding of the general features of the bounce solution from what we learned in 
Case A. 
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First, the bubble at the moment of nucleation would be a Q-Ball of the min- 
imum size, where the surface energy is as important as the volume energy. The 
charge density p* inside the Q-ball would be much larger than the initial charge 
density, po- The minimum size of a Q-ball could be rather large if w* is very close 
to m so that the energy gain by the charge is small and so a lot of charge is needed 
to compensate the surface energy. We know that there will be a large current 
will flow into inside the bounce wall from outside in this case because the interior 
charge density pi = ffv% is much larger than pq. With a similar argment given 
to Case A, the current would be large inside and outside wall at the poles, and 
would be zero outside and large inside the equator. The energy density outside the 
bounce wall at the pole is larger than that outside the bounce wall at the equator. 
The barrier energy at the poles would in turn be lower than that at the equator, 
and so the wall tension at pole will be lower than that at the equator. Thus, the 
bounce solution in Case C would also be elongated along the r direction in Case 
A. Fig. 6 shows such a bounce solution for Q-ball nucleation. 

Once a Q-ball is nucleated, it will grow but very slowly. The reason is that a 
Q-ball can grow only when it swallows the charge from outside and that there is 
not much charge around it because the formation itself have already diluted the 
initial charge density around its neighborhood. This is consistent with the picture 
that the bounce solution is elongated along the r direction, which also implies a 
slow terminal velocity as argued after Eq.(4.9). The explicit nature of the Q-ball 
nucleation and expansion would however require the better analysis and would not 
be attempted here. 
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5. Discussion 



We have studied the phase transitions in the theory of a complex scalar field 
with a global U(l) symmetry when there is non zero initial charge density We 
have discussed the metastability condition on the possible initial configurations 
and the possible inhomogeneous final configurations. We argued that there are 
many cases of the phase transitions to be studied in the theory. We have set the 
Euclidean formalism of the bubble nucleation when there is nonzero charge density. 
We applied our formalism to the case when the initial charge density is small and 
when the phase transitions involve the symmetric phase as the initial configuration 
or a part of the finial configuration. Here we studied the characteristics of the 
bounce solutions and the bubble evolution. Our system is shown to have a rich 
variety of the possible phase transitions, and could be a good simple toy model of 
the phase transition involving charges, the supersonic bubble wall, and the sound 
wave. 

However, there are still many loose ends and questions we have not attempt 
here. One of the interesting questions seems what is the later development of the 
bubbles. Depending on the cases of the phase transitions, there is a possibility 
of rich dynamics. Additional interesting questions to be explored are about the 
phase transition between the asymmetric phases and about the phase transition 
when the initial charge density is not small. Finally, we note that it is rather 
straightforward to extend our formalism in Sec. 3 to the case involving the nonzero 
local gauge charges. 
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FIGURE CAPTIONS 



1) Plot of e(p, /) and U(f) for a typical case. 

2) Plot of two potentials Ui(f) and U 2 (f). 

3) Plot of y/2U(f)/f 2 for three potentials U h U 2 , U 3 . 

4) The bounce solution for Case A. The dashed circle is the wall of the 0(4) 
symmetric bounce solution. The dotted lines indicate the charge flow. The 
solid ellipse is the wall of the deformed bounce. 

5) The bounce solution for Case B. The dashed circle is the wall of the 0(4) 
symmetric bounce solution. The dotted lines indicate the charge flow. The 
solid ellipse is the wall of the deformed bounce. 

6) The bounce solution for Case C. The dotted lines indicate the charge flow. 
The solid ellipse is the wall of the bounce for Q-ball nucleation. 
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